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ABSTRACT 

We present Space Telescope Imaging Spectrograph (STIS) spectral images 
of the HH 30 stellar jet taken through a wide slit over two epochs. The jet is 
unresolved spectrally, so the observations produce emission-line images for each 
line in the spectrum. This rich dataset shows how physical conditions in the jet 
vary with distance and time, produces precise proper motions of knots within 
the jet, resolves the jet width close to the star, and gives a spectrum of the 
reflected light from the disk over a large wavelength range at several positions. 
We introduce a new method for analyzing a set of line ratios based on minimizing 
a quadratic form between models and data. The method generates images of the 
density, temperature and ionization fraction computed using all the possible line 
ratios appropriately weighted. In HH 30, the density declines with distance from 
the source in a manner consistent with an expanding flow, and is larger by a 
factor of two along the axis of the jet than it is at the periphery. Ionization in 
the jet ranges from ~ 5% — 40%, and high ionization/excitation knots form at 
about 100 AU from the star and propagate outward with the flow. These high- 
excitation knots are not accompanied by corresponding increases in the density, 
so if formed by velocity variations the knots must have a strong internal magnetic 
pressure to smooth out density increases while lengthening recombination times. 
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INTRODUCTION 



Included with the first catalog of Herbig-Haro objects (IHerbig I Il972l ). HH 30 is in 
many ways the prototype for all jets from young stars. HH 30 is located n ear the L1551 
molecular cloud in the Tauru s star formation region at a distance of 140 pc (IKenyon et al. 



1994 ). iMundt &: Friedl (119831 ) first noticed the jet-like nature of HH 30 from ground-based 



images, and identified both a main jet and a counterjet. The system achieved a great 
deal of fame when images from the Hubble Space Telescope (HST) revealed a well-resolved, 
flared reflection nebu la on either side of an opaque circumstellar disk oriented nearly edge-on 
( jBurrows et al.lll996l . hereafter B96), and a collimated clumpy jet emerging exactly perpen- 
dicular to the disk plane on both sides of the disk. The HST images made it clear that 
flared disks do not collimate high-velocity optical jets on large scales; instead, jets become 
collimated within a few tens of AU of the disk plane, leading to the modern inter pretation 
of magnetically-dominated jets driven by accretion disks (e.g., 



Ferreira et al.l 120061 ) 



The reflection nebulae at the top and bottom of the flared disk in HH 30 are bright, 
and quickly became the primary testbeds for models of scattered light within accretin g 
protostellar disks and envelopes (IWood et al.l Il998l ; ICotera et al.l l200ll ; IWood et al.ll2002l ). 
The models show the disk inclination to be 82 degrees, with the brighter (north) side of the 
flow tilted towards the Earth. The nebulae vary significantly in brightness and in morphology, 
reflecting conditions present at the base of the accretion disk where material fa lls onto the 
star fB96. IWood fc Whitnev!ll998l : IWood et~aI1l200ol : Iwatson fc Stapelfeldtbooj ). Like other 
stellar jets, HH 30 extends to parsec scales from the source (ILopez et al.lll996l ). The jet is not 
straight, but curves, possibly as a result of precession (B96). Seeing-limited ground-based 
observations typically do not resolve individua l knots in HH 30 w ell. Jet widths inferred 
from deconvolved ground-based images (~ OA "; IMundt et al.lll99lT) are about twice as large 
as those found closer to the star with HST by lRay et al. Jl996h . while B96 report the jet to 
be unresolved everywhere even with HST. 

Because the HH 30 flow is in the plane of the sky and produces a reliable source of 
new emission-line knots every few years, it is an ideal target for studying jet physics. The 
nearly edge-on disk blocks the stellar light effectively, so one can in principle trace the jet 
to within a few tens of AU of the source and determine how physical conditions such as the 
density, temperature, ionization fraction and collimation behave as the jet first emerges from 
the accretion disk. The orientation of the flow close to the plane of sky means there are no 
significan t projection effects in im ages. More importantly, the internal velocity dispersion 
is small (lAppenzeller et al.ll2005l ). and there are no strong bow shocks in the flow, so the 
line widths in the jet are unresolved with typical grating spectrometers. Hence, a spectrum 
taken through a slit wide enough to include the entire jet gives images of the jet in all the 
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emission lines with a single exposure. Moreover, one can remove the continuum light by 
subtracting background at wavelengths adjacent to each emission line in the spectrum. 

In this paper we present a set of deep (10 — 30 ksec) images of the HH 30 jet in all the 
emission lines visible in the jet from 3500A through 1/zm. The new data are dithered, and 
make it possible to follow the collimation in the jet down to ~ 20 AU of the source. The 
images were taken at two epochs, so we can also measure precise proper motions for the first 
time in this flow. We present the first high-quality spectrum of the reflection nebula from 
the near-UV through the near-IR, and we quantify how colors vary with distance from the 
source in the reflection nebula. 

But the main focus of this paper is on the physical conditions in the jet, which we can 
only explore in detail by having many line ratios to use as diagnostics. Determining the best 
solution for the physical parameters of a plasma from an observed set of non-independent line 
ratios is a common problem in nebular astrophysics, and we discuss an analysis technique 
that includes all the information present in the spectrum. The resulting images of the 
ionization fraction, density, and temperature of the jet at each position give new insights 
into how jets are heated as they emerge from the source. 

2. OBSERVATIONS AND DATA ANALYSIS 

Our program used the Space Telescope Imaging Spectrograph (STIS) on HST to observe 
HH 30 in the fall of 2000 (hereafter epoch 1) and again in the fall of 2002 (hereafter epoch 2; 
Table 1). The low-resolution gratings G430L and G750L were both used in epoch 1, and we 
reobserved HH 30 with G430L in epoch 2. The medium resolution grating G750M was used in 
both epochs. The HH 30 reflection nebula was considerably brighter in the first epoch than in 
the second epoch. The second epoch G430L observations are most useful for quantifying how 
the reflection nebula changed with time, while the two epochs with the G750M grating are 
ideal for measuring proper motions in the resolved jet knots, which emit brightly in the red 
lines of [O I], [N II], [S II] and Ha. The G430L observations were taken over a time interval 
of about three weeks in the fall of 2000, but the source showed no significant variability over 
this period so we combined these observations into a single exposure. Likewise, the G750M 
frames taken 14 Aug and 22 Aug 2002 are essentially identical, and were coadded. 

Standard pipeline processing is inadequate for these data because it leaves too many hot 
pixels in the images and loses spatial resolution during the rectification process. We reduced 
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the data using tasks in the C ALSTIS package of STSDAS within IRAF. @ A customized IRAF 



IMFORT task described by lHartigan et al.l (120041 ) removed hot pixels. Flux-calibrated 2D 



images have units of surface brightness (erg cm _2 s _1 A _1 arcsec -2 ), and we checked to make 
sure that our flux calibrations matched those produced by pipeline reductions. As a check of 
the flux calibration, the surface brightness in epoch 1 of the brightest HH knot in the [O I] 
6300A G750L spectrum agreed to better than 10% with the G750M observations. 

Each set of observations for a given epoch and grating consists of at least four sep- 
arate exposures dithered by fractional pixels to improve the spatial resolution. The driz- 
zle task within the STS DAS package of IRAF combined individual images into composites 



( IFruchter Sz Hookil2002l ). The plate scale in the final composites is 0.02536 arcseconds per 
pixel, or 3.55 AU per pixel for a distance to HH 30 of 140 pc. The spatial resolution of 
HST is close to the diffraction limit at the wavelengths of the bright red emission lines in 
HH 30 (~ 0.07 arcseconds at 6500A). After dithering, the G750L spectra have a dispersion 
of 2.441A per pixel and cover 5234A - 10230A, the G750M spectra have 0.277A per pixel 
between 6290A and 6857A, and the G430L spectra have 1.374A per pixel from 2907A to 
5728A. The signal-to-noise-ratio (S/N) in the reflection nebula is low blueward of ~ 3200A 
and redward of ~ 1/zm, so we truncated the spectra accordingly. 

The slit was 52 arcseconds in length and aligned at PA —147.35° to lie along the HH 30 
jet. With the exception of one set of images taken through the 0.2 arcsecond slit to check for 
internal velocity dispersion (see below), all data were taken through the 2 arcsecond wide slit. 
We performed a blind offset from a nearby star to center HH 30 in the slit. The closest star 
suitable for this offset was over 2 arcminutes away, which could introduce pointing errors of 
a few tenths of an arcsecond. However, these errors are small compared with the size of the 
slit, and observations of the scattered light in the disk show that the slit was centered on the 
source to within ~ 0.2 arcseconds. The C-shaped reflection nebula that surrounds HH 30 
extends about 2.5 arcseconds perpendicular to the slit, and therefore fits almost entirely 
within the slit. Because the nebula also reflects emission lines, an image of the reflection 
nebula appears near the base of the jet for each of the strong emission lines. Figure 1 shows 
a portion of the G750M spectrum. 

Observations of jets taken through slits that are wider than the jet produce emission-line 
images at each line, provided the internal velocity dispersion of the emitting gas is unresolved. 
In our case the slit is ~ 5 times wider than the jet (Fig. 1). Internal velocity widths within 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. The Space Telescope Science Data Analysis System (STSDAS) is distributed by the 
Space Telescope Science Institute. 
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HH jets are nearly always < 100 kms 1 , and are unresolved with the low resolution gratings 
G430L and G750L. However, we require the medium-resolution grating G750M to image 
lines separated by 10 - 20 A, such as [S II]A6716 and [S II]A6731, and [N II]A6583 and 
Ha. With most jets, wide slit observations with G750M would be difficult or impossible to 
interpret because spa tial and spectral i nformation are mixed. However, HH 30 lies nearly in 
the plane of the sky (IWood et al.ll2000l ). and has no strong bow shocks near the source that 
splatter material alon g the line of sight. Hence , emission-line velocities and widths in the jet 
near the star are low. lAppenzeller et al.l (120051 ) found that while scattered light generated a 
weak emission tail that extends to ~ 100 kms -1 in the spectrum of the disk of HH 30, this 
emission is down in flux by a factor of ~ 100 from that of the jet , where the emission -line 
profile lies entirely within a 40 kms -1 interval. Coude spectra by iMundt et al.l (119901 ) also 
show narrow line widths in the jet close to the star. An emission-line point source will be 
imaged over ~ 2.5 pixels in the dithered 750M observations, corresponding to ~ 0.7A, or 
~ 40 kms -1 , so the jet knots should be unresolved in our spectra. To check this, we also 
imaged the jet with the 0.2 arcsecond slit, and found the emission lines to be unresolved. 
As a further check, the G750L and G750M emission-line images of well-separated lines such 
as [O I] 6300 and [O I] 6363 are identical, which demonstrates that the jet is unresolved 
spectrally in both gratings. 

To create an image of a particular line, we simply extract a box centered on that line. It 
is straightforward to measure the continuum from the reflection nebula on either side of the 
emission line and subtract it from the image (e.g. lHartigan et al.l 120041 ) . Images in various 



emission lines appear in Figure 2. The relative positions in the x-direction (perpendicular 
to the jet) are in principle fixed by the wavelength solution. However, the scatter in that 
solution of ± 0.3 pixels is improved by a factor of ~ 2 by shifting the x-positions of the 
images to align the brightest knot in the jet. The jet knot used for alignment appears fairly 
symmetrical in all the lines, so the uncertainty in its x-position is ~ 0.1 pixel. 

Aligning each emission-line image in the y-direction (along the jet) is more involved 
because the star in HH 30 is not visible, no other stars lie along the slit, and emission lines 
have offsets from one another in the y-direction owing to heating and cooling along the flow. 
Fortunately, HH 30 has a bright reflection neb ula visible on either side of the disk, which 
appears as a dark lane in images of the source (IWatson Stapelfeldtl 120041 ). and it is easy 
to measure the position of this lane in our spectral images. 

At our request, K. Wood and M. O'Sullivan generated a series of ID slices of a model 
of the HH 30 reflected light at several wavelengths between 3000A and 9500A. The model 
slices are perpendicular to the disk through a 2-arcsecond slit, and show that the location 
of the minimum of the reflection nebula at all wavelengths remains fixed to within a half 
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pixel relative to the projected location of the star. The models show that the star is located 
0.092 ± 0.010 arcseconds away from the minimum of the reflected light. We measured the 
y-location of the minimum of the reflection nebula adjacent to each emission line, and shifted 
each image accordingly. 

Spectra extracted from the flux-calibrated surface brightness images or ratios between 
sections of such images must be corrected for the point-spread function (PSF) of STIS, which 
worsens at the reddest wavelengths. The calibration is accomplished with a photometric 
correction table. The values affect the relative line ratios by < 10% for wavelengths < 
8500A, which includes all of our lines except for the [C I] 9850 doublet, where the correction 
is ~ 40%. 



3. RESULTS 



3.1. Reflection Nebula 



The reflection nebula at the base of the HH 30 jet has been known for many years and 
has been the subject of several investigations. Most recent studies mak e use of HST obser- 
vations at optical (B96. IWatson & Stapelfeldtl 12004 ) and near-infrared (jCotera et al.ll200ll ) 
wavelengths to understand the scattering properties and spatial distribution of the dust. The 
surface brightness of the nebula varies by as much as 1.5 magnitudes at optical wavelengths, 
and there is evidence that reflected light fro m hot spots on the stellar surface causes most 
of th e observed asymmetries in the nebula (IWood & Whitney! 1 19981 ; IWatson fe Stapelfeldt 
20041 ) . Further support for this idea comes from observations of em ission lines in the nebula , 
which have broad velocity tails that are likely to be scattered light (lAppenzeller et al.ll2005l ). 



The reflection nebula is easiest to study in the low-resolution (G430L and G750L) data 
because the signal-to-noise in the continuum is higher than it is for the medium-resolution 
data. A single emission line reflected from the disk appears in the slitless data as an image 
of the disk in that line, and occurs at some level in all the bright lines in Figure 2, but is 
especially bright in the Balmer lines, which come in part from the star. Each row of the 
slitless images that crosses the disk (the jet lies along the columns) provides a spectrum of 
the reflected light that falls along that row, although all absorption and emission lines will 
appear broadened by the physical size of the disk in the extracted spectrum. One can obtain 
an integrated spectrum over a portion of the disk by simply summing the appropriate range 
of rows. 



Spectra of the entire disk show that the total flux of the reflected continuum varies by 
over a magnitude at V, consistent with previously published results. The blue spectra from 
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13 August 2000 and 3 September 2000 (Table 1) are nearly identical, so we combined them 
into a single epoch. The nebula became fainter in the second epoch of blue data (18 August 
2002) by a factor that increases steadily from about 2.0 at 3500A, to 3.0 at 5500A (Fig. 3). 
A strong Balmer jump and many emission lines are present in the blue spectra. The red 
spectrum taken 25 October 2000 was intermediate in flux between the two blue spectra in 
the region where the wavelengths overlap. 

The emission-line images in Fig. 2 show that the broad base to the emission-line profiles 
in Fig. 3 come from the disk, while the narrow cores arise in the jet. The ratio of the broad 
to narrow components, measured from the medium-resolution red spectra (the lines are not 
blended there), is about 0.9 for all the forbidden lines and ~ 4 for Ha. Both the forbidden 
lines and Balmer lines from the disk are likely to be scattered light, with the star being the 
primary source of illumination for the Balmer lines and the jet the primary source for the 
forbidden lines. 

The Ha image of the disk in Fig. 4 shows an asymmetry that appears in both epochs of 
the G750M data, with the Ha brighter on the left (northwest). The asymmetry persists in 
the G750L images, and so is not caused by spectrally resolved scattered blueshifted emission 
from the outflow. The asymmetry is not present in the forbidden li nes, indicating that it 



origin ates in the strong Ha emitting region in the vicinity of the star. IWatson fe Stapelfeldt 



( 120041 ) noted similar asymmetries and attributed them to hot spots on the stellar photosphere 
caused by accretion columns from the circumstellar disk. The Ha scattered light images 
appear truncated on the right (southeast), but this position likely marks the edge of the slit, 
because this truncation appears in the forbidden lines as well. 

In addition to the Balmer jump and emission lines, there are weak TiO absorptio n 



bands at 7800A and 7100A in the disk spectrum (Fig. 2, see also lAppenzeller et al.l 120051 ) 



but the signal-to-noise is too low to constrain the spectral t ype of the star well. Usin g a 



STIS spectrum of the M3.5 weak-line T Tauri star HBC 358 (IHartigan fc Kenyonll2003l ) as 
a photospheric template, the veiling is about 4 at these wavelengths, while the K7 wTTS 
Lk Ca7 (primary) gives a veiling of 1.5. By chance, the brightness of the nebula in epoch 
1 agrees to within a few tenths of a magnitude with the broadband V and I magnitudes 
reported by B96. 



Watson fc Stapelfeldtl (120041 ) found that HH 30's reflection nebula became bluer further 
out in the disk, and we also see this effect in our data. The ratio of the G430L spectrum 
of the outer half (~ 1/3 arcsecond) of the blueshifted reflection nebula to a spectrum of the 
inner half (closest to the star) rises in a nearly linear fashion from ~ 1.0 at 5500A to ~ 1.3 
at 3500A. 
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3.2. Proper Motions 

Knots in the blueshifted (northeastern) half of the HH 30 jet exhibited substantial 
proper motions of ~ 15 pixels in the 1.8 years between epoch 1 and epoch 2. To measure 
these motions we combined the continuum-free images of the five brightest forbidden lines 
visible in the G750M spectra - [O I] 6300, [O I] 6363, [N II] 6583, [S II] 6716 and [S II] 6731 



into a single image for both epochs. We then used the code described in lHartigan et al. 



(120011 ) to measure proper motions. The results appear in Table 2 and in Figure 5. Velocities 



in the blueshifted jet are rather low for an HH jet, ~ 130 kms 1 . Velocit y differences of 



30 kms 1 are present, similar to those observed in other HH flows (e.g., lHartigan et al. 



200lh . 



Velocities of knots in the redshifted counterjet are higher than they are in the main, 
blueshifted jet. Three knots are present in the first epoch at distances of 0.63, 0.96, and 2.38 
arcseconds from the star, which we label D, E, and F, respectively. Knots E and F both have 
proper motions of ~ 225 kms -1 . The motion for knot D is less certain because knots exist 
in epoch 2 at distances 1.04 and 0.68 arcseconds from the star. Assuming the outer knot is 
knot D, the proper motion is ~ 150 kms -1 . It is possible that knot D had not emerged fully 
from the high extinction region of the disk plane in epoch 1, which would result in a lower 
proper motion for this object. In any event, the redshifted portion of the flow is ~ 1.7 times 
faster than its blueshifted counterpart. Having o ne side of a flow much fas ter than the other 



side is not unusual for stellar jets (e.g., RW Aur; lLopez-Martin et al.ll2003l ). However, in the 
case of HH 30 it does not appear that the knots were ejected at the same time (Table 2). 
There are knots in the HH 30 blueshifted jet at larger distances from the source than knot 
C in Fig. 5, but these are either partially or completely excluded by our slit because the jet 
curves at larger distances. The jet appears slightly tilted by about 2 degrees with respect 
to the position of the slit, and the motions of the knots all lie along the axis of the jet to 
within the uncertainties of the observations. 

It is instructive to compare our results with previous proper-motion measurements. B96 
determined proper motions in HH 30 with WFPC2 from broadband red images separated by 
about a year. While HST can easily detect proper motions over this time interval, the B96 
measurements have large errors because the first epoch exposure in particular was not very 
deep, and the images include bright continuum from the disk which makes it more difficult 
to measure accurate positions for knots within an arcsecond of the star. Proper motions of 
the six knots in the blueshifted jet measured by B96 have a large scatter, from as low as 54 
kms -1 to as high as 258 kms -1 , with an average of 142 kms -1 . Ground-b ased images do 



not resolve individual knots well close to the source, but iMundt et al.l (J1990j) measured 170 



± 50 kms 1 for the proper motion of a distant knot ~ 10 arcsesconds from the star. 
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The combination of deep exposures (10 — 30 ksec), a factor of four smaller pixel scale 
in the dithered STIS images compared with WFPC2, precise continuum subtraction, and 
improved analysis software greatly reduces the errors in the proper motions to ± 4 kms -1 . 
With this precision it is possible to distinguish real velocity differences between adjacent 
knots and to predict where the features we observe now should have appeared in the January 
1995 images of B96 (the spatial resolution is too low in the ground-based images to make a 
useful comparison). The knot labeled 'A' in the blueshifted jet in Fig. 5 was ejected after 
1995, but knot 'B' should have been located at a distance of 0.40 ± 0.03 arcseconds from the 
star in January of 1995. This position corresponds reasonably well, given the uncertainty 
of the source position and the different filter used, with knot 9501-N, which B96 reported 
at a distance of 0.51 arcseconds. However, the B96 tangential velocity of only 54 kms -1 is 
much lower than our value of 149 ± 4 kms -1 (Table 2). The difference may be caused by 
uncertainties in measuring proper motions within the reflected light cavity of the disk with 
the broadband B96 images. 

Our object 'C, which we measure as having a proper motion of 125 ± 4 kms -1 , should 
have been located 2.11 ± 0.03 arcseconds from the star in January 1995. B96 report three 
knots in this 9502-N (1.1 arcseconds, 258 kms -1 ), 9503-N (2.0 arcseconds, 158 kms -1 ), 
and 9504-N (2.7 arcseconds, 84 kms -1 ). It appears that these three objects visible in 1995 
have now merged into a single structure - object C extends for ~ 1.8 arcseconds, about the 
distance between 9502-N and 9504-N in 1995. 

Some structural changes are evident in Fig. 5. Knot A became more elongated between 
2000 and 2002, while knot B widened and faded, and the portion of knot C closest to the 
star also faded somewhat. These changes are likely caused by internal motions, shocks, 
rarefactions and cooling still unresolved in these images. However, internal motions in these 
knots cannot exceed ~ 40 kms -1 to remain consistent with observations of narrow intrinsic 
linewidths in these knots. A new knot (labeled 'N' in Figs. 6 and 7) was emerging from the 
source during our observations, but proper-motion measurements are not possible because 
the feature has not yet detached clearly from the source. 

3.3. Jet Collimation 

The factor of four finer spatial scale of our dithered slitless STIS images relative to 
narrowband WFPC2 images, the lack of continuum, better hot pixel removal strategies, and 
long exposure times create an unprecedented database for measuring the width of the HH 30 
jet. The best images for this measurement are [S II] 6716 and [S II] 6731 with G750M, and 
[O I] 6300 with G750M and G750L, because in each of these cases the lines are bright and 
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unblended, and there is no stellar emission-line component. We measured the FWHM by 
extracting spatial profiles every three pixels along the jet in each emission-line image and 
fitting a Gaussian to the profile with the SPLOT command in IRAF. Errors are dominated 
by uncertainties in the baseline, and are larger near the star where the disk reflects the 
line emission in the jet. The instrumental FWHM for slitless observations is that of the 
PSF of STI S. We used the exper imentally measured values and uncertainties of the PSF 
reported in lHartigan et al.l (120041 ) to deconvolve the FWHM, and propagated the errors of 
the measured FWHM with the uncertainties in the deconvolution kernel to obtain errors for 
each deconvolved point. 

Results are shown in Fig. 6 for epoch 1, which have both G750L and G750M data, 
and in Fig. 7 for epoch 2. The widths measured from the G750L and G750M exposures in 
[O I] 6300 are identical, so the jet is unresolved spectrally in all images. The upturn in the 
linewidth for the [O I] points < 15 AU is likely to be instrumental in origin. The brightness 
of the jet drops sharply inside 20 AU, implying less contrast with the reflected emission lines 
in the disk. Also, there will be some light in the wings of the PSF from the bright jet at > 
30 AU that contaminates these measurements. 

The collimation properties of the HH 30 jet have been the subject of some debate in 
the literature. Using HST, B96 observed an opening angle of about 3 degrees, and measured 
the width of the jet to within about 70 AU of the star. They concluded that th e jet was 
unresolved at this distance, with a FWHM < 20 AU. In contrast. iRay et al.l (119961 ) analyzed 
the same image and found the jet to be spatially resolved everywhere, with a FWHM ~ 
35 AU at a distance of 50 AU from the star, implying a very wide opening angle for the 



j et at the source. Deconvolving ground-based images with an FFT algorithm, iMundt et al. 
( 119911 ) reported an even larger FWHM of ~ 60 AU at a distance of ~ 350 AU, although this 



measurement is very difficult to do because the jet is much narrower than the seeing disk 
(which is ~ 1 arcsecond or 140 AU FWHM). 

Our new observations settle this controversy. Fig. 6 shows that the jet is indeed spatially 
resolved e v erywh ere by a t least 5cr, though the width is only about half that reported by 
Ray et al.l (119961 ) and by IMundt et al.l (Il99ll ). and is near the upper limit of B96. The jet 
widens gradually from 14 ± 3 AU at 20 AU from the source to 36 ± 4 AU at 500 AU, implying 
an opening half-angle of 2.6 ± 0.4 degrees, consistent with the opening angle determined by 
B96 at larger distances. 



The HH 30 jet is similar to other jets (e.g. lHartigan et al.l |2004j) in that it does not 
project to a point at the source. However, this observation does not mean that the jet 
emerges from a disk of radius 14 AU. Rather it is more likely that the opening angle of the 
flow is wider near the source, as expected for a wind launched from an accretion disk. There 
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is no obvious correlation betwe en the presence of a bright knot and the width of the jet, as 
also noted by iRay et al.l (119961 ). However, there is some structure in the width of the jet, 
including relatively wide areas at about 230 AU and 310 AU during epoch 1 (Fig. 6). These 
areas propagated downstream about 50 AU by epoch 2 (Fig. 7), at about the jet velocity. 
Apparently these sort of irregularities in the jet develop near the source and the flow simply 
carries them along (see also Fig. 5). 

On the fainter, redshifted side we coadded all the forbidden lines together to improve 
the signal-to-noise but we were only able to measure jet widths within ~ 400 AU of the star 
(Fig. 8). The measurements are too uncertain to constrain the opening angle well, but the 
jet is clearly resolved everywhere and is wider than it is on the blueshifted sid e. At larger 
distances from the star, previous images of the region (IMundt et al.lll990l . Il99ll . ; B96) have 



also shown that the redshifted flow is wider and has a larger opening angle than its blueshifted 
counterpart. Being able to compare the opening angles, velocities, and widths in both a jet 
and a counterjet is a potentially powerful diagnostic tool for analyzing the physics of how 
jets are launched provided enough sources can be observed to derive meaningful statistics. 

For a freely expanding flow, the observed opening half-angle corresponds to a Mach 
number of 21.8. Using a jet velocity of 130 kins -1 we find a sound speed of 6.0 kms -1 , 
which corresponds to a temperature of only 2600 K for a ratio of specific heats 7 = 5/3 and 
mean molecular weight 1. This temperature is far too low to explain the fluxes in bright lines 
that originate from upper levels such as [S II] 4068. The more likely value is ~ 10 kms" 1 , 
which should occur for low-excitation forbidden lines in a mostly neutral cooling zone of a 
shock. 

The observed opening angle is consistent with the temperature if there is a source of 
confining pressure that is on the order of the thermal pressure. The pressure cannot be 
provided by an ambient external medium, because shear between the jet and this medium 
would rapidly heat the interface between the two fluids, creating strong shock waves in each 
that do not appear in the observations. Instead, the most likely candidate for a confining 
pressure is a toroidal magnetic field. Equating the magnetic pressure with the thermal 
pressure, using a total density of 10 6 cm -3 , and sound speed of 10 kms -1 (T = 7260 K) we 
find the magnetic field B ~ 5 mG at ~ 300 AU, i n agreement wi t h the 'fiducial' values of 
field strengths in jets at this distance discussed by lHartigan et al.l (120071 ). 
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3.4. Emission Line Ratios within Individual Jet Knots 

Fig. 9 shows the observed ratio images of the four bright, red forbidden lines for the 
two epochs. The signal-to-noise is excellent in the ratio images, with uncertainties in the 
ratios of < 5% except within ~ 50 AU of the source at the edges of the jet, where reflected 
light introduces higher uncertainties. The effect of scattered light is highest on the blue 
(left) side of the ratio images that involve [N II] 6583, because scattered Ha from the disk 
contaminates that side of the emission-line image. The images were registered to the stellar 
position (marked with a horizontal white line) as described in section 3.1. The shorter 
horizontal lines in the figure mark locations of bright knots N, A, B, and C in the summed 
forbidden line images for each epoch, and are included as guides. 

The easiest emission-line ratio to interpret is that of [S II] 6716/ [S II] 6731, which 
decreases monotonically with increasing electron density. This ratio image clea rly shows that 



the electron density increases close to the source (see also lBacciotti et al.lll999l ). although the 
ratio is in the high density limit close to the star and therefore only gives a lower limit to the 
electron density in those areas. The [S II] ratio image also shows that the density declines 
near the edges of the jet. This decline has been inferred from STIS mapping in other jets 



(.e.g. 



Bacciotti et al.ll2000i ). but Fig. 9 shows the effect with unprecedented spatial resolution 
and clarity. The intensities in the combined forbidden line images in the two left panels 
are far more concentrated along axis than is the electron density. Two effects contribute to 
this behavior: first, the volume emissivity is proportional to n 2 so small density fluctuations 
translate into larger flux variations; and second, the brightness in the emission-line images 
integrates along the line of sight, making the jet appear brighter along its axis even for a 
uniform density cylinder. 

The second set of image ratios in Fig. 9 is [N II] 6583/ [S II] 6716+6731. This ratio 
depends primarily on the ionization fraction, but also has a weak positive dependence on 
both the electron density and temperature. Hence, this ratio follows 'excitation' such as 
shocks, increasing wherever the gas becomes hotter, denser, and more ionized. The N II/S II 
ratio images from the two epochs are fascinating, and show substructure not apparent in 
the individual images. The ratio is high in the two bright knots A and B, and these high 
excitation regions move outward with the flow velocity. The highest excitation portions 
of the jet lie closer to the source than the brightest portions of the knots in both cases. 
Additional structure in this ratio appears in the second epoch in the vicinity of knot N. 

There is a region of relatively low [N II] /[S II] near the source, also present in the 
[N II]A6583/[0 I]A6300 ratio (not shown). This decrease appears to be real and not instru- 
mental. A registration error along the jet of 2 pixels between the two images would produce 
a similar falloff of [N II] /[S II] near the source, but the reflected light profiles in the original 
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images show that registration uncertainties are at most 0.5 pixels. As a final check, the 
[N I]A5199+5201/[N II] A6583 ratio also rises < 30 AU from the source, indicating that a 
drop in ionization fractio n occurs there (see nex t section). Beyond knot B, the [N II] /[S II] 



ratio declines as noted by iBacciotti et al.l (119991 ) 



The final ratio shown in Fig. 9 is [O I]A6300/[S II] A6716+6731, which is primarily sensi- 
tive to density, although the [O I]A6300 line also has a positive dependence on temperature 
not present in the [S II] lines because the collision strength of [O I] is a function of tem- 
perature. In a mostly neutral plasma like a stellar jet, most oxygen is O I, and essentially 
all sulfur is S II. The upper levels for the [O I]A6300 and [S II]A6716+6731 transitions have 
similar energies so the temperature dependence of the ratio is weak. However, the critical 
density of [O I]A6300 (~ 10 6 cm" 3 ) is significantly higher than that for [S II] A6716+6731 
(~ 10 4 cm -3 ), so [O I]A6300 becomes stronger relative to [S II] A6716+6731 close to the 
star where the density is high. Comparing the [S II]A6716/[S II] A6731 ratio image with 
[O I]A6300/[S II]A6716+6731 shows that the latter ratio increases sharply as the S II lines 
approach the high density limit, so together these three emission lines constrain the density 
well everywhere in the jet. 



4. Discussion 

4.1. Physical Conditions Along the HH 30 Jet 

A main motivation of this work is to convert images of emission-line ratios like those 
in Fig. 9 to images of T e , N e and X# = Nhii/O^hii+^hi) at each epoch. In this section 
we describe how to create these images, consider how these variables change along and 
perpendicular to the jet, and investigate the sensitivity of results to assumptions such as 
abundance values and charge exchange coupling. The results give the first clear picture of 
how jets are heated as they emerge from their sources. Our data can also be used to infer 
the reddening and mass loss rates but we defer discussion of these issues to a second paper, 
where we also include blended blue doublets such as [S II]A4068+4076 and [O II]A3727+3729 
into the analysis. 



4-1.1. Determination of T e , N e , and Xh 

We use images of the four brightest forbidden lines obtained with the G750M setting, 
[O I]A6300, [N II]A6583, [S II]A6716, and [S II]A6731, to find T e , N e , and X H at each point 
in the jet where the S/N in the ratio between any pair of these lines is > 5. N I, N II, O I, 
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O II, and S II each have five levels populated in nebular conditions: 3 P, 1 D, 1 S for I and 
N II and 4 S, 2 D, 2 P for II, N I, and S II. We do not use [0 I]A6363 or [N II]A6548 for 
anything other than flux calibration. Ratios of [0 I]A6363 / [0 I]A6300 and [N II]A6548 / 
[N II] A6583 are constant to within the errors of measurement everywhere along the jet, and 
equal the ratios of their respective A-values, as expected from lines that originate from a 
common upper level. The other bright red line in the HH 30 jet is Ha, but this line has a 
prominent reflected stellar component that contaminates the jet emission. 

We constructed a model level population for N I, N II, 01,0 II, and S II for a grid of N e 
and T e with the atomic parameters summarized in Appendix B. With solar-like abundances 
(0=8.82, N=7.96, S=7.30 on a log scale with H=12.0) we calculated emission-line ratios for 
\h = — 1 in increments of 0.01, log N e = 2 — 7 in steps of 0.1, and T e /10 4 K = 5 — 
25 in steps of 0.5 assuming the ionization of N and O are tied to H by charge exchange as 
described below, and all S is S II. The latter assumption makes sense as no [S III] or [S I] 
lines are present in the spectrum. The best fit physical parameters were defined at each 
point by minimizing the quadratic form C described in Appendix A over this grid of 209100 
models (100x51x41). 

The resulting images of N e , T e and X# appear for each epoch in Fig. 10, where, as in 
Fig. 9, the horizontal white line marks the location of the star and short white lines denote the 
position of the knots in the jet at each epoch. Dividing the electron density by X# produces 
a total density image, and the product of T e and X# gives an 'excitation' image that is 
sensitive to hot ionized gas, as one should find behind a shock wave (Fig. 11). The density 
images largely resemble what appears in the [S II] image ratio in Fig. 9, with the dominant 
behavior the decline of density with distance from the source, as h as been observed before 



at lower spatial resolution in the 1-D plots of iBacciotti et al.l (119991 ). Density measurements 



are still possible close to the source where the red [S II] lines are in the high density limit 
because the best fit solution also includes line ratios that involve [O I]A6300, which is not in 
the high density limit. 

The decline of the total density with distance shown in Fig. 12 resembles that produced 
by a radial conical flow that originates from a finite source region. The jet is also a factor 
of two denser along the axis than it is at its edges, and the bright portions of the jet are 
typically the densest. The high-density knots clearly move along the jet, and the higher 
S/N images of the second epoch show what appear to be density filaments at the limit of 
the spatial resolution in knot B. When interpreting these images one should keep in mind 
that the observations integrate along the line of sight, which tends to smooth over density 
gradients in the jet. 

The images of X# and excitation show a rapid rise, followed by a gradual decline with 
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distance from the source (see also lBacciotti et al.lll999l ). Figs. 11 and 12 show that this varia- 
tion is not smooth; instead, the high excitation regions show sharp, almost linear boundaries, 
and move outward with the flow. The regions of high excitation occur on the side of the 
knot closest to the source for both knots A and B. In knot A we see that this area of high 
excitation appeared between the first and second epochs, indicating a heating event. 

The rise of the ionization/excitation in the first 100 AU as the jet emerges from the 
source in a series of well-defined heating events is a critical observational result of this paper 
that affects all models of jet launching, so it is important to verify that the result does 
not depend upon the emission-line ratios used in the analysis. As noted above, using other 
lines close to the source typically requires reddening corrections, but in this case we can 
use the G430L [N I]A5199+5200 image to confirm the low value of X# found close to the 
source from the red line ratios. If the ionization fraction remained constant, the [N II] A6583 
/ [N I] A5 199+5201 ratio should increase steadily towards the source because the critical 
density of [N II] A6583 exceeds that of [N I]A5199+5201, and because the reddening increases 
toward the source. However, images of this ratio in the last panel of Fig. 11 clearly show 
that the ratio is lower close to the source, which can only happen if the ionization fraction 
is lower there. 



4-1-2. Effects of Charge Exchange Coupling, Reddening, and Abundances 

The assumption that charge exchange sets the ionization fractions of O and N given 
the ionization fraction of H is worth investigating further, as it is required to solve for 
the hydrogen ionization fraction and temperature from a database that consists of only 
[N II]A6583, [O I]A6300, [S II]A6716, and [S II]A6731. The charge exchange cross section 
between H and O is large and has almost no temperature dependence because the ionization 
energies of H and O are nearly the same. However, the charge exchange cross section is 
three orders of magnitude lower for H-N than it is for H-O, so there is some question as to 
how well the ionization fractions of N are tied to those of H, es p eciall y since the ionization 



potentials of H and N differ by ~ 0.9eV. iBacciotti & Eisloffell (119991 ) argue that the H-N 



charge exchange rate is high enough to fix the (H I/H II) /(N I/N II) ratio as long as the 
hydrogen ionization fraction is < 50%, but their model did not include photoionization, 
which affects H and N differently owing to the different ionization potentials. 

To test the validity of H-N charge exchange coupling we created a time-dependent 
photoionization model, varying the ionizing spectrum and strength, and included the atomic 
physics of collisional ionization, recombination, dielectronic recombination, photoionization 
and charge exchange (see Appendix B for references) to determine how closely the ionization 
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fraction ratios of H and N follow those predicted by charge exchange coupling. 

The first case we tested was to singly ionize all the H, N, and O, and follow the ion- 
ization fractions as the gas recombined. As expected, the (H I/H II)/(0 I/O II) ratio 
approached rapidly to the value predicted from charge exchange (to within 1 part in 10 4 
after 0.1 trec(H), where trec(H) is the H recombination timescale), and remained equal to 
the charge exchange ratio as the gas recombined. The ionization fraction of N was 10% lower 
than the charge exchange value after 0.5 trec{H) improving to 5% lower after 2 trec{H)- 
Experiments with input blackbody photoionizing spectra yield similar results as long as the 
ionization fraction of H is < 90%. Even the extreme case where the photoionizing flux is low 
and H starts out completely ionized with N completely neutral (so that charge exchange is 
solely responsible for ionizing N) equilibrates to within ~ 20% of the charge exchange pre- 
diction after one trec{H). We conclude that the H-N charge exchange coupling assumption 
is well-justified. 

Another issue is whether or not we should attempt to deredden the fluxes of the four 
emission lines. In principle one might use the Balmer decrement (Hce/H/3) for this purpose. 
The observed Balmer decrement in the bright knots varies between about 3.0 and 3.5, and 
rises to > 5 in the reflection nebula and between the bright knots. This variation of the 
Balmer decrement is likely to be caused by differing decrements in the scattered light from 
the star and the intrinsic emission from the jet. Hence, we cannot use the decrement di- 
rectly to find the reddening, and instead must incorporate various blue doublets such as 
[O II]A3727+3729, [S II]A4068+4076, and [N I]A5199+5201 into the analysis. We defer this 
work to a subsequent paper. The Balmer decrement in the bright jet knots indicates that 
A v < 1, so differential reddening between [O I]A6300 and [S II]A6731 is a factor of < 1.04, 
and between [N II] A6583 and [S II] A6731 a factor of < 1.02, small compared with the un- 
certainties in the relative abundances of O, N, and S in the Taurus star forming region. We 
ignore these small effects of differential reddening for the current analysis. 

The density maps are largely independent of the abundances because the [S II]A6716/[S II] A6731 
ratio strongly constrains N e when log N e < 4. For a fixed set of observed emission-line ratios, 
reducing the nitrogen abundance by a factor of two increases X# by a factor of 1.8 and T e 
by a factor of 1.2. Overall, the fit to the line ratios becomes worse for about 3/4 of the pixels 
and better for the remaining 1/4. Hence, it is possible to constrain abundances with these 
data but the uncertainties are relatively high. We also defer this analysis to a subsequent 
paper. 
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5. Implications for Jet Heating 

The images presented in Figs. 9—11 provide the best information to date concerning 
how jets are heated as they emerge from their sources. The observations show that some 
process heats and partially ionizes the jet from < 10% within ~ 20 AU of the star to as high 
as 30% -40% at ~ 100 AU from the star. The ionization increase is not a steady process, 
but instead produces relatively sharp boundaries, and areas of high ionization move outward 
with the flow velocity. The heating is not caused by photoionization from the source, because 
the ionization fraction is initially low close to the source, and increases outward. Similarly, 
a stationary focusing shock in the jet does not explain why the jet produces distinct knots 
with relatively high ionization fractions, and there is no indication of any narrowing in the 
jet that should accompany such a shock. 

The most obvious candidates for heating and partially ionizing the flow are shock waves 
produced by supersonic variations in the flow velocity. This scenario explains why the jets 
are heated at some distance from the source, and internal shocks produce sharp ionization 
features that propagate with the flow velocity, in agreement with the observations. However, 
one important aspect of the data does not fit this model — the density should increase 
substantially in the postshock regions where the excitation is high, and it does not. 

It may be possible to produce sharp rises in the ionization fraction without increasing 
the density substantially if the shocks are magnetic. In this case the postshock pressure is 
dominated by the magnetic field, and relaxes on the timescale for Alfven wave propagation, 
while the ionization in the shocked region relaxes on the recombination timescale. Non- 
steady MHD jet models with atomic physics accurate enough to follow ionization fractions 
have not yet been done, but the ionization signature of a supersonic magnetic collision should 
linger substantially longer than the density enhancement does in stellar jets. Magnetic field 
strengths in jets at distances of ~ 100 AU from the source are poorly understood, but recent 
models of time-d ependent MHD pulse s indicate that the Alfven velocity could be > 50 km/s 



in these regions (lHartigan et al.ll2007l ). so an Alfven wave would cross the jet on a timescale 
of < 2 years. In comparison, the recombination time for H gas with a density of 10 4 cm -3 is 
longer, ~ 10 years, approximately equal to the flow time from knot A to knot C. 

The general decline of ionization with distance means that more heating events occur 
close to the star (as expected for stochastic velocity perturbations). Beyond ~ 1000 AU 
from the star, proper-motion measurements and high-resolution HST images of jets (e.g. 



Hartigan et all l200ll ; iHeathcote et al] Il996l ) show that jets are heated by internal shock 
waves. At these distances, small velocity perturbations have washed out, so the jet gradually 
recombines and remains mostly neutral and cool until it encounters a major working surface. 
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6. SUMMARY 

We have constructed deep optical emission-line images of the HH 30 stellar jet using 
STIS spectroscopy taken through a wide slit. The system has a fortunate geometry, with 
the axis of the jet nearly in the plane of the sky. Hence, it is possible to observe the jet 
without projection effects and we can trace the flow to within ~ 20 AU of the star. Low 
radial-velocity dispersions within the flow mean that the STIS spectra are effectively slitless, 
and produce images for each emission line, including those for which HST has no narrowband 
filters. The resulting data set, dithered for maximum spatial resolution and observed at two 
epochs, enables us to diagnose physical conditions in the region where the jet first emerges 
from the accretion disk. 

Individual knots in the HH 30 jet show distinct proper motions between the 1.8 years 
that separate the two epochs. Motions in the blueshifted jet range from 116 kms -1 to 149 
kms -1 , with uncertainties of ± 4 kms -1 . The jet has a resolved spatial width of FWHM ~ 
14 AU at a distance of 20 AU from the source, and has a constant opening half-angle of 2.6 
degrees. The narrow opening angle is consistent with the observed velocity and sound speed 
if there is an additional confining magnetic pressure from toroidal fields in the jet. Velocities 
in the redshifted counterjet are ~ 100 kms" 1 higher than they are in the main jet, and the 
counterjet is less well-collimated than the main jet. 

Spectra of the reflected light from the HH 30 disk reveal a veiled, late-type photosphere 
with both permitted and forbidden emission lines and a prominent Balmer emission jump. 
The reflected light from the disk varies substantially in both morphology and in brightness, 
and is bluer at larger distances from the star. 

We developed a new analysis technique to find the best fit for the electron temperature 
T e , electron density N e , and hydrogen ionization fraction X# given a set of observed and 
model emission-line ratios. The method, based on minimizing a quadratic form, has the 
advantage that by construction it uses all the information present in each available line 
ratio, and appropriately weights the fit of each ratio by the uncertainty. 

Our analysis focused on six ratios involving the four bright nebular lines of [N II] A6583, 
[O I]A6300, [S IIJA6716, and [S II]A6731, with additional ratios using [O I]A6363, [N II]A6548, 
Ha, H/3, and [N I]A5199+5200 used to verify the results. The images from each epoch produce 
the first high-resolution images of T e , N e and X# in a stellar jet. The density in the jet is 
highest close to the source, and declines in a manner similar to that of radial conical flow 
from a finite source region. The density in the jet is larger by a factor of two along the axis 
than at its edges. Distinct bright knots are also denser, and propagate at the flow velocity. 

Maps of the ionization fraction and excitation (defined as the product of T e and X H ) 
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show that the jet emerges with an ionization fraction of only ~ 10% , which increases to ~ 
30% at ~ 100 AU from the source. Regions of locally higher ionization propagate with the 
flow, and exhibit sharp spatial boundaries. At least one such knot formed between the two 
epochs. Surprisingly, the high-excitation knots are not accompanied by a density enhance- 
ment, suggesting that they may originate from velocity variability in a highly magnetized 
flow. 

The two epochs reported in this paper were an accidental benefit of a delay in HST 
scheduling. Stringent ORIENT requirements limited visibility windows for the spacecraft, 
and NICMOS took priority for all usable 2001 dates. As a result, our data were taken in two 
narrow time intervals, making it possible to study the motions and variability of the jet and 
the disk. Stellar jets are dynamical systems, and we must watch them move to be sure we 
understand the physics that governs them. Because of the two epochs, we now know where 
the jet is heated and we are beginning to understand the process that causes the heating. 
Multiple epoch observations with the highest spatial resolution possible are clearly the best 
way to study these remarkable objects. 
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7. APPENDIX A: Measuring Electron Densities, Electron Temperatures and 
Hydrogen Ionization Fractions From Emission Line Ratios 

In this appendix we describe our method for extracting electron densities, tempera- 
tures, and hydrogen ionization fractions from a set of emission-line ratios and uncertainties. 
In general, an observed emission-line ratio from an optically thin plasma where the domi- 
nant processes are collisional excitation, collisional deexcitation and radiative decay depends 
on the reddening, the relative abundances of the two elements, their respective ionization 
fractions, the electron temperature and the electron density. One can use the flux ratio 
of lines with closely-spaced upper energy levels such as [S II]A6716/[S II] A6731 to measure 
the electron density directly, because these ratios are independent of abundances, redden- 
ing, and ionization fractions, and insensitive to the electron te mperature. This method has 
been applied to HH objects for at least a half century now (jBohml Il956l ). Other line ra- 
tios that arise from the same ion but where the upper levels differ markedly in energy (e.g. 
[O I]A5577/[0 I]A6300]) depend on both the temperature and density, so with the density in 
hand one can measure the temperature (IBrugel et al.lll98ll ; lHartigan et al.l 120041 1 . 



A variant of the above ideas introduced by iBacciotti fc Eisloffell ( 119991 ) was to use the 
additional constraint of strong charge exchange coupling between H and O, and H and N, 
which ties the H I/H II ratio to O I/O II and to N I/N II given the temperature of the 
gas. Assuming standard abundance ratios of N/O, O/S and N/S, and taking all the S to be 
S II, one can then use [S II]A6716/[S II] A6731 to obtain the electron density as before, and 
employ [N II]A6583/[0 I]A6300, [S II]A6716+6731/[0 I]A6300, and [O I]A6300/[N II]A6583 to 
estimate the temperature and ionization fraction. Because these lines all have similar upper 
level energies, the temperature determination is the most uncertain, and relies mostly upon 
the temperature dependence of the collision strength of [O I] for the measurement. The 
main advantage of assuming that charge exchange ties the ionization fraction of hydrogen 
to that of nitrogen and oxygen is that one can then estimate the desired physical quantities 
from bright lines in the less-extincted red part of the spectrum that are close enough in 
wavelength to make differential reddening corrections negligible. 

Estimates of N e , T e , and X# made in this manner implicitly assume that these parame- 
ters do not vary substantially over a size scale equal to t he spatial resolution of the observa- 



tions. However, as shown clearly by IBrugel et al.l (ll98ll ). and more recently by iPodio et al. 



(120061 ). temperatures and densities of HH objects measured with different line ratios pro- 
duce markedly different results. The reason for this behavior is that HH objects are shock 
waves. Emission lines from shocks radiate at different places in the cooling zones, with 
high-excitation lines and lower critical densities more prominent near the front, and low- 
excitation, high critical density lines stronger at larger distances. Unless the cooling zones 
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are resolved spatially, almost never the case with ground based observations of HH jets, 
constructing line ratios when the component lines sample different physical conditions can 
lead to inconsistent results. Hence, great care must be taken when interpreting large scale 
trends in HH shocks when the coolin g zones are unresolved . Fortunately, HST resolves the 



cooling zones for most knots in jets (IHeathcote et al.l 119961 ). so it is possible to apply the 
above analysis to many such images. 

There are several aspects to the above methods of analysis that are less than ideal. The 
[S II] A6716/[S II] A6731 ratio is usually in the high density limit close to the source, and 
provides only a lower limit to the electron density there. Density information is contained in 
other line ratios, but the analysis does not use all the ratios. Moreover, it is unclear how to 
incorporate additional line fluxes into the analysis. For example, if in addition to the four 



l ines (7S II]A6716, [S II]A6731, [N II]A6583, and [O I]A6300) used bv Bacciotti fc Eisloffel 



( 119991 ) we also observe [N I]A5199+5201, what is the best way to include this observation in 
the analysis? In principle, the new observation defines four new emission-line ratios, each 
with a different uncertainty. 



7.1. Definition of the Problem 

Given a set of n observed line fluxes Fi,F 2 ,...F n , what is the best way to determine the 
physical parameters N e , T e , and Xh = n(H II)/n(H I+H II)? 

We begin by discarding the flux normalization. Absolute line flux measurements from 
ground-based spectra are rather rare, as they require photometric conditions and are com- 
plicated by atmospheric dispersion and pointing uncertainties. Physically, the volume of 
emitting gas along the line of sight and the reddening determine the absolute fluxes for these 
optically thin lines. In HH jets, which are inherently clumpy, the absolute fluxes determine 
the volume filling factor of the emitting material if one knows the width of the jet, but the 
absolute fluxes do not affect estimates of T e , N e , and X#, which depend only upon the flux 
ratios. Hence, our analysis deals only with flux ratios. The extra step of determining filling 
factors is easy to add after one finds T e , N e , and X#, assuming the spectra are flux-calibrated. 

Most studies present line ratios relative to some normalization, typically H/3 or Ha = 
100. However, there is no reason a-priori to favor any particular line for the normalization. 
Consider Table 3, which shows artificial data for two models and an observation of three 
emission lines, Ha, [O I]A6300 and [N II] A6583. At first glance it appears that model B fits 
the observations better than model A because the [N II] flux is off by a factor of two in both 
models, but model B gets the [O I] right while model A is again low by a factor of two. 
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However, there is an additional line ratio to consider, which is [0 I]/[N II], and model A 
predicts this ratio correctly while model B is off by a factor of two. Hence, if we consider all 
line ratios equally, both model A and model B are equally good fits, because they both get 
one ratio correct and miss the other ratios by factors of two. 



Motivated by the above discussion, we seek a measure that will use all the emission-line 
observations in equal, and appropriately weighted manners to determine the best fit set of 
the parameters N e , T e , and X#. Given n line flux measurements we can construct p = 
n(n— 1)/2 distinct pairs. The deviation of the model from the data should be the same for 
a ratio and for its inverse, and to ensure this behavior we define a new set of p ratios r^, 
where r^ = ln(Fj/Fj) for some integers i, j = 1, 2, n, with i < j. The most obvious choice 
to determine a best fit model is to minimize the quadratic form 



where m k is the model prediction of the line ratio k from a given model, and a k is the 
observational uncertainty in the line ratio. 

The quadratic form C has all the properties we desire for a line ratio fitting algorithm. 
Each of the emission-line ratios contributes on an equal basis to the value of C according to 
the uncertainty in the ratio. The best fit is determined from line ratios and not fluxes, and 
information from each ratio is automatically incorporated into the fit for the best physical 
parameters. If, for example, ratio r» is unaffected by the temperature then it will contribute 
equally to the value of C for all models that differ only in the temperature. Line ratios that 
give only upper limits will do so in the values of C as well. Adding new emission lines into 
the analysis is trivial. Applied to models A and B in Table 3, we see that Ca = C_b, as 
desired. We define the 'best' fit to a given set of parameters T e , N e , and X H to be the one 
that minimizes C for a set of observed line ratios. 



7.2. The Quadratic Form C 




(Al) 



7.3. The Probability Distribution of C 



An outstanding issue is how one should interpret the values of C. Although this exercise 
has no bearing on the choice of the best fit, which by definition is the one with the lowest 
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value of C, it is still useful to investigate the probability distribution of C because in a formal 
sense the values of C determine confidence intervals for rejecting a null hypothesis that the 
model exactly matches the data. 

For astronomical spectra it is reasonable to take the individual r& to be distributed 
normally. Equation Al looks a lot like a Xn distribution, but it is not because the r& are not 
independent of one another. For example, in Table 3 any two line ratios suffice to determine 
the third. If m and r are the p-dimensional vectors of a set of model ratios and observations, 
respectively, and the r^ have a normal distribution with mean nij and variance of, then a 
standard result from multivariate statistics is that the probability density function 



(2*) 



p/2 



exp 



mFV- 1 ! 



m ) 



(A2) 



is a multivariate normal, where V is the covariance matrix of r, and the superscript 'T' 
denotes a transpose. U sing the above notat ion, the quadratic form H 



r— m 



V 



111 ) 



is distributed as a x\ (IHoog & Craig! Il995l ). In the case of independent r i; the covariance 
matrix is diagonal and we recover the standard definition of x\- While we could also minimize 
H to find the best fit, with the benefit that H then has a known probability distribution, 
doing so requires inverting the covariance matrix at each point within an emission-line map, 
an extra computational step that did not seem warranted given that the quadratic form C 
exhibits all the properties we desire for a best fit solution to the problem. 

The probability distribution of C = x T Ax has been studied tho roughly in the math- 
ematical literature. Examples include independe nt normal variables (iCochranl I1934J ) . nor- 
mal variables with nonzero cova riance (IBoxi 119541) . and arbitrarily distributed vectors with 
nonzero covariance when A = I (IBlacher 1 120031 ) . When, as in our case, the Xj are normally- 
distributed with nonzero covariance, the probability distribution of C = YJ% j=i-^ij x i x j * s 
given by a weighted sum of x 2 distributions, each with 1 degree of freedom: 



P(C) = AiX? (A3) 

where the weights Aj are the eigenvalues of VA, with V the covariance matrix. For our 
application, Ay = Sija^ 2 . Hence, the probability distribution differs for each set of line 
ratios depending on the uncertainties in the ratios of the emission lines at that point. For 
the simplest case of three line ratios with equal uncertainties o and equal covariances Oy, 
setting Ay = 5ij<7~ 2 , b = o^/o 2 and keeping terms of order b (0 < b < 1) in a power series 
expansion, we find the three eigenvalues Aj = 1, l±\/36. The eigenvalues approach unity 
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and the distribution becomes xl as b^O, as expected for independent variables. 

Experimenting with typical line ratios and uncertainties, the probability distribution 
of C for p-ratios typically resembles a Xp f° r large values of C, but the actual distribution 
depends on the input data. For example, were we to include a very low S/N emission line 
with a large uncertainty into a list that initially has four emission lines, the four new line 
ratios would not significantly increase the value of C even though p increases from 6 to 10. 
The best fit model would not change, as desired, because the new noisy line does not affect 
C significantly. 



8. APPENDIX B: Atomic Parameters for Emission Line Models 



Energy levels for N I, N II, O I, O II, and S II derive from the compilation of IMendoza 



(119831 ) . We used published values of the Einstei n- A coefficients for the forbidden transi 
tions within N I (butler & Zeippenl Il98j) , N II (bell et alJ Esaj; Istorey. & Zeippenl 5)00 



Froese-Fischer fc Saha 



(IMendoza fc Zeippenl Il982h . 



r (z Zieippenllib)84i). in ii imeii et aun^a ; storey &; z,erpi 
1983), O I (jBaluja &: Zeippenl Il988l ). O II (jZeippenl Il982l ) . 



and S II 



W e interpolated collision strengths for N I transitions from tabular values of lBerrington fc Burke 
Jl98lh . where we fit transitions to the 4 S level with log 10 (fi) = a + b*logi (T/10 4 ) + 



c*logio(T/10 4 ) 2 , with a, b, and c con stants. Internal fine structure ratios in N I were 
taken from Table II of iPradhanl (11976 ). Col lision strengths within the 2 P and 2 D levels 
used the power law fits of jDopita et al.l (119761 ). Collision stren gths for the N II lines com e 
from iHudson fc Belli (120051 ) . Collision strengths for O I are from iBerrington fc Burkd (119811 ) , 



except for transitio ns within the 3 P levels, which arise from power law fits we made to the 
tabulated va lues o f iLe Dourneuf &: Nesbetl (119761). In a rec ent study of planetary nebulae, 



Wang et a. 



of 



Zeippen 



(2004) found the O II A- values of IZeippenl (119821) to b e more accurate than those 



( 1987). a nd th e collision strengths of 



Pradhanl (119761 ) to be better than those of 



McLaughlin &: Belli (119931 ). Collision strengths for the lines of S II are from iKeenan et al. 
dl996h . 



Radiative recombination coeffi cients are from iPeauinot et al.l (1199 ll ). and dielectronic 



recombination coefficients a r e from iNussbaumer fc Storey! (119831 ) at low temperatures and 
Landini &: Monsignori Fossil (119901) at h igh temperatures. We used the charge-excha nge rates 
compi led by lKingdon fc Ferland ( 1996 ). the collis ional ionizati o n rate coefficients of Voronov 



( 119971 ). and the photoionization cross sections of IVerner et al.l (119961 ) 
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Table 


I. STIS Observing Log 


Grating 


Number of 


Total Exposure 


Dates 




Exposures" 


Time (sec) 




G430L 


5 


13323 


13 Aug & 3 Sept 2000 


G430L 


4 


10538 


18 Aug 2002 


G750M 6 


1 


2146 


3 Sep 2000 


G750M 


4 


10467 


25 Oct 2000 


G750M 


12 


31404 


14 Aug & 22 Aug 2002 


G750L 


4 


10262 


25 Oct 2000 



a Each exposure was CR-split, and used the 52x2 slit except as 
noted. 

b 52x0.2 slit. 



Table 2. Proper Motions in the HH 30 Jet 



Object" 


Separation 6 


A X c 


A Y c 


Velocity^ 


Age 6 


A 


0.78 


1(2) 


-175 (5) 


116 (4) 


4.5 


B 


1.70 


-1 (3) 


-224 (5) 


149 (4) 


7.6 


C 


3.20 


-7 (2) 


-188 (5) 


125 (4) 


17.0 


D 


-0.63 


19 (5) 


-227 (5) 


151 (5) 


2.8 


E 


-0.96 


31 (10) 


-325 (10) 


217 (9) 


3.0 


F 


-3.20 


13 (9) 


-343 (10) 


229 (9) 


6.9 



a Object boundaries are defined in Figure 5. 

b Distance in arcseconds from the center of the box that de- 
fines the object to the HH 30 star for epoch 1. A negative 
number corresponds to the redshifted counterjet. 

c Proper motions (and errors) in milli-arcseconds per year. 
The X and Y-axes are defined as in Figure 5, with the negative 
Y-axis corresponding to motion along the jet for the blueshifted 
portion of the flow. 

d Velocity in the plane of the sky in kms" 1 for a distance of 
140 pc. 

e Age of the knot in years relative to 2000.82 UT, the first 
epoch of observations. 



Table 3. Example of Observed and Model Fluxes for Thr ee Lines 
Line Observation Model A Model B 



Ha 
[O I] 6300 
[N II] 6583 



1.00 (0.10) 
0.40 (0.04) 
0.20 (0.02) 



1.00 
0.20 
0.10 



1.00 
0.40 
0.40 




H H 30 1 35 AU = 0.25" = 1 pixels 




Fig. 1. — Spectral images of the HH 30 system with the G750M grating. The locations 
of the star, disk plane, and bright emission lines are marked. Reflected light is present at 
the top and the bottom of the opaque disk, which is nearly edge-on in this object. The 
brightest (blueshifted) portion of the jet extends below the disk (to the northeast) in this 
and subsequent figures, while a fainter redshifted counterjet exists above the disk plane (to 
the southwest). The slit is much wider than the jet, so these data generate an image for 
each emission line. 




Fig. 2. — Epoch 1 emission- line images of the HH 30 jet. The slit was oriented at a po- 
sition angle of -147 degrees. Images of [O II]A3727, [S II]A4070, H/3, [N I]A5200, and 
[O I]A5577 are from G430L, while images of [O I]A6300, Ha, [N H]A6583, [S II]A6716, and 
[S II]A6731 are from G750M, and images of [O II]A7325 and [C I]A9850 are from G750L. 
Adjacent contours are spaced by a factor of two, with the lowest contour level (in units of 
10 -i5 ergcm-VWsec- 2 ) equal to 3.47 for [O II]A3727 and [S II]A4070, 2.17 for H/3 and 
[O I]A5577, 1.73 for [N I]A5200, 1.93 for [O II]A7325, 3.08 for [C I]A9850, and 0.70 for all the 
other lines. Line identifications labeled with an asterisk are doublets or quartets that are 
partially resolved spectrally. The horizontal white line marks the location of the source. 
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Fig. 3. — Spectra of the HH 30 Disk. Top: A medium-resolution spectrum taken 14 Aug 
2002, obtained by integrating the entire light of the disk and any emission from the jet 
that is superposed atop the disk. The broad base to the profiles results from emission lines 
reflected from the disk, and not from high-velocity material. The narrow emission comes 
from spatially confined knots in the jet. Bottom: Low-resolution spectra of the entire HH 
30 disk. Vertical offsets between the spectra are caused by flux variability between the 
observations. Weak TiO bands are present in the red spectrum around 7100A and 7800A. 




Fig. 4. — An Ha image of the HH 30 disk and jet taken with the G750M grating on 25 Oct, 
2000. The disk appears in scattered light, and has a bright extension to the left. The slit 
partially truncates the disk on the right side. Adjacent contours are separated by a factor 
of y/2, and the lowest contour is 2.2 x 10~ 15 erg cm~ 2 s~ 1 arcsec~ 2 . 




Fig. 5. — Proper motions in the blueshifted component of the HH 30 jet. The objects labeled 
A, B, and C move ~ 0.4 arcseconds (56 AU) over a time interval of 1.8 years. Contours 
are spaced by a factor of in these [S II] images, with the lowest contour 1.1 x 10~ 15 erg 
cm~ 2 s _1 arcsec~ 2 . The source position is denoted with an 'X'. A new knot is also emerging 
from the source. 
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Fig. 6. — Top: Measured FWHM of the jet as a function of distance from the source in four 
emission line images taken during epoch 1. One arcsecond corresponds to 140 AU. The two 
horizontal lines mark the upper and lower limits to the measured instrumental FWHM. Each 
point represents an average of three lines. Bottom: Deconvolved jet widths. The jet width 
gradually increases with distance from the source. The curve is the relative intensity of the 
forbidden lines with distance, obtained by adding together all the bright red emission-line 
images and coadding the fluxes for each row. Knots labeled A, B, and C are used to measure 
proper motions. A new knot 'N' is emerging from the source. Extinction from the nearly 
edge-on disk obscures the jet at distances < 20 AU from the source. There is no correlation 
of the FWHM with the presence or absence of a bright emission knot. 
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Fig. 7. — Same as Fig. 6 but for epoch 2. Irregularities in the jet width have propagated 
outward with the jet velocity since the first epoch in Fig. 6. 
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Fig. 8. — Same as Fig. 6 but for the redshifted counterjet using the sum of [0 I] 6300, [O I] 
6363, [N II] 6583, [S II] 6716 and [S II] 6731. The Gaussian FWHM for each point was 
measured from a spatial cut perpendicular to the jet, coadding five rows (0.125 arcseconds) 
for each point. The jet is wider on the redshifted side than it is on the blueshifted side. 
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Fig. 9. — Emission line images and ratios for epochs 1 and 2. Horizontal white lines mark 
the locations of the bright knots N, A, B, and C. Each of the three emission-line ratios has 
its own scale bar, color-coded for that ratio. The [S II] image ratio shows the jet is denser 
along its axis and closer to the source. Areas of high [N II] /[S II] and [O I]/[S II] on the 
side closest to the source in knots A and B move outward with the jet, and both ratios show 
a gradual decline with distance from the source, whose position is denoted by a horizontal 
white line. 



[0 I] 6300+6363 + [N 111 65S3 Electron Density Electron Temperature H Ionization Fraction 

+ [S II] 6716+5731 " 

epoch 1 epoch 2 epoch 1 epoch 2 epoch 1 epoch 2 epoch 1 epoch 2 




Fig. 10. — Same as Fig. 9 but for the electron temperature, electron density and hydrogen 
ionization fraction X# = H II/(H I+H II). The area of high excitation (high X# and T e , 
color coded as yellow and white) in knot B moves outward with the flow. The jet emerges 
with a very low (< 0.1) ionization fraction, and becomes more ionized beyond ~ 50 AU. At 
large distances the ionization declines. 
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Fig. 11. — Same as Figure 9 but for the total density, excitation, and [N II] 6583/[N I] 5200 
ratio. The total density rises to a value of ~ 4 x 10 5 near the source, and is higher (red 
and yellow colors) in the bright knots and along the axis of the jet. Dense filaments, most 
easily seen in the epoch 2 image, are present within knot B. The excitation image, created 
by the product of the ionization fraction and the electron temperature, shows several areas 
of heating between ~ 50 AU and 500 AU from the star. The [N II] /[N I] ratio, which is not 
dereddened, declines close to the source even though the density and reddening both increase 
toward the source, and both factors act to increase the ratio. The only explanation for the 
decline of the ratio close to the source is that the jet is more neutral there, in agreement 
with the ionization map in Figure 10. 
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Fig. 12. — Density structure in the HH 30 jet. Top: The electron density and the total 
density decline with distance from the source. The mean and standard deviation of the 
electron density and total density in a box that extends 3 pixels along the axis of the jet and 
9 pixels perpendicular to the jet produce each point and its errorbar. Density fluctuations 
within the flow move along at the jet velocity between epochs. The solid curve depicts a 
conical flow emerging from a base with a radius of 10 AU. Bottom: Density perpendicular to 
the flow axis in epoch 2 calculated over three distance ranges from the source. The density, 
calculated from the ratio of the N e and X# images, is larger by about a factor of two on axis 
relative to the edges of the jet. 



